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Abstract 

The quantum electromagnetic dielectric tensor for a multi species plasma is re-derived from the 
gauge invariant Wigner-Maxwell system and presented under a form very similar to the classical 
one. The resulting expression is then applied to a quantum kinetic theory of the electromagnetic 
filamentation instability. Comparison is made with the quantum fluid theory including a Bohm 
pressure term, and with the cold classical plasma result. A number of analytical expressions are 
derived for the cutoff wave vector, the largest growth rate and the most unstable wave vector. 

PACS numbers: 52.35.Qz, 41.75.-i, 03.65.-w 
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I. INTRODUCTION 



In the last years there has been an increasing interest on electromagnetic (as opposed 
to purely electrostatic) quantum plasma phenomena. With the emergence of new fields 
like spintronics |l| where magnetic effects are fundamental, it is essential to have a deeper 
knowledge on electromagnetic quantum plasma models. For instance, recently Eliasson and 
Shukla |2| have reported on the Bernstein modes in a magnetized, degenerated quantum 
plasma described by the Vlasov- Maxwell system. In this work, the fermionic nature of the 
charge carriers was taken into account thanks to the underlying Fermi-Dirac type equilib- 
rium. Consequently, the upper hybrid dispersion relation for a degenerate quantum plasma 
was derived. However, the employed kinetic equation was classical (Vlasov), so that the 
second quantum effect besides statistics, or quantum diffraction effects, was not included. 
To take care of the typical quantum phenomena like wave packet spreading and tunneling, 
one must resort to a quantum kinetic equation for the reduced one-particle Wigner function, 
which is the quantum analog of the usual classical reduced one-particle distribution func- 
tion in phase space. For these reasons, it is advisable to go one step further and formulate 
a general kinetic electromagnetic theory for linear waves in quantum plasmas. There are, 
nevertheless, already many instances where linear and nonlinear wave propagation in elec- 
tromagnetic quantum plasmas were treated. For instance, one can cite the quantum Weibel 
instability, the dense plasma magnetization by the electromagnetic waves, the temporal dy- 
namics of spins in magnetized plasmas, stimulated scattering quantum instabilities and the 
analysis of self-trapped electromagnetic waves in a quantum hole, as reviewed in Ref. jsl by 
Shukla and Eliasson. 

In the present work we show a kinetic treatment of the quantum filamentation instability, 
significantly generalizing the earlier fiuid-based model {4, 5|. Moreover, the basic kinetic 
model we use is the evolution equation for the gauge-invariant Wigner function [6;], which 
fully takes into account quantum diffraction. In addition, we assume a counter-streaming 
zero-temperature Fermi-Dirac equilibrium, so that the quantum statistical effects due to the 
exclusion principle are also incorporated. 

This work is organized as follows. In Section II, we derive the general electromagnetic 
dispersion relation for linear waves in quantum plasmas. In Section III, the quantum kinetic 
filamentation instability is worked out. The dispersion relation is compared with the results 
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from the classical zero-temperature fluid model and the quantum hydrodynamic model, with 
a Bohm potential term. Also the behaviors of the cuttof wave vector, of the most unsta- 
ble wave vector and the corresponding largest growth rate are studied. Further, analytic 
estimates for these quantities are provided. Section IV is reserved to the conclusions. 



II. DIELECTRIC TENSOR 



The electromagnetic dielectric tensor of a coUisionless classical plasma can be derived in 
a ,u,te sta.da.d n..nne, f.on, the Vlasov-Maxwel. system QB. In a ve,v sioHa. wa. the 

quantum version of the same tensor can be derived from the Wigner-Vlasov system, where 
the Wigner equation \^ is the quantum equivalent to the Vlasov one. Following Ref. lo| |. 
the gauge invariant Wigner function /j(r, v,t) for particles of species j with charge qj and 
mass rrij, obeys a Vlasov-like equation introduced by Stratonovich {g]. This same equation 
was put into an illuminating form by Serimaa [11], according to 
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where B = B(r, t) and E = E(r,t) are the magnetic and electric flelds respectively. To 
compute Avj, E^ and B^, one has flrst to Taylor-expand in powers of h the electromagnetic 
fields in the integrands and then perform the integrals. For instance, up to second-order in 
h, one has 
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The kinetic equation ([T]) follows from the von Neumann equation solved by the reduced one- 
body density matrix. As it stands, Eq. ([1]) is manifestly gauge-invariant since it depends 
only on the fields. 



Assuming an equilibrium Wigner function /j'(v) in a zero equilibrium electromagnetic 
field, one can set 



/. = /» + //(v) exp[z(k-r-c^t)] 
E = El exp[i(k - r — ut)], 
B = Bi exp[i(k ■ r — ut)] , 



(7) 
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where fj{'v),'Ei and Bi are first order quantities. In this case it follows from Eqs. (jSHH]) 
that 



Ej — E Lj , Bj — B Lj , 



defined in terms of the operators 
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The operator Lj in expression ( ITTj) is understood in a Taylor-expanded sense. For instance 
up to second-order in h we have 



L, = l + ^(k-^)2 + ... . 



(12) 



Linearizing Eq. ([T]) we get 
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while the linearized Maxwell- Faraday and Ampere-Maxwell equations resp. gives 



k X El = wBi , 

k X Bi = -ifio ^ qjUoj j d\wf] - ^Ei 



(14) 
(15) 



where /xq is the vacuum permeability constant, c is the speed of light and riQj is the equihb- 
rium number density of specie j. The normalization j dv = 1 is applied. The line of the 
calculation is then very similar to the classical case [12]. Eliminating and Bi between 
Eqs. (HSD-dlS]) and using Eq. ^ yields T(k, w) ■ Ei = 0, with 
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where Upj = {nQjq'j / {mjeo)y^'^ is the plasma frequency for specie j. 

The expression (JT6l) can be simphfied using the two following properties, 
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which can be proven by series expansion of Lj besides an integration by parts. In this way 
we find 



T(k, u) = -^e{k, w) + k ® k - r I 



where 
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From Eq. ( IT9|) the general dispersion relation for linear electromagnetic waves in quan- 
tum plasmas can be written as det(TQ,^) = 0. As is well known, linear dispersion relations 
for quantum plasmas have been discussed for decades. For instance, Lindhard [l3| obtained 
both the longitudinal and the transverse dielectric tensors for a quantum electron gas in a 
Fermi-Dirac equilibrium. For arbitrary equilibria, the transverse and longitudinal dispersion 
relations for quantum plasmas have been found by Klimontovich and Silin jl4| and Bohm and 
Pines 15 1 respectively. In addition Silin and Rukhadze 16 1 and Kuzelev and Rukhadze |l3] 
have found the general dielectric tensor in Eq. (!20!) in a slightly different presentation, us- 
ing the gauge-variant Wigner-Maxwell and Schrodinger-Maxwell systems respectively. Also 



Kelly [18( has found the dispersive properties of a quantum plasma under an homogeneous 
magnetic field, using a gauge-dependent Wigner formalism. Furthermore, the straightfor- 
ward derivation based on the Stratonovich-Serimaa gauge invariant Wigner-Maxwell system 
should be compared to the more cumbersome calculations needed in the gauge-dependent 



formalism 



18|, applicable only in the particular case of homogeneous magnetic fields. In 



addition, Eqs. (fT9l |20|) are valid for multi-species quantum plasmas and arbitrary equilibria, 
encompassing both longitudinal and transverse perturbations in an unified way. 
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FIG. 1: (Color Online) Distribution function considered for two identical counter streaming beams 
with here VJ, = 2Vf- 

III. THE FILAMENTATION INSTABILITY 
A. Dispersion Equation 



We now turn to investigate the filamentation instability with wave vector k = [k, 0, 0), for 
a beam-plasma system where the flow is aligned with the y axis. For two symmetric identical 
counter streaming beams at T = with densities nb/2 and velocity Vb, the dispersion 
equation reads [19 1 

Eyy = k'^ C^/U'. (21) 

As can be seen in Figure [H the distribution function F is the sum of two spheres of radius 
Vp (Fermi velocity) centered around ±Vf,y, 

3 



(22) 



and = outside the Fermi sphere. For this equilibrium a substitution of variables readily 
simplifies Eyy given by Eq. fl2U]) into, 



2U2 



-yy 



dv 



47r Vl^u' J^^vp (ou - kv^f - h'ky{4m' 
This quadrature can be calculated switching to spherical coordinate with 

vcose, {v,e) G [0,Vf] X [0,7r] 
f COS*/? sin 6*, G [0, 27r] 
w sin sin 6'. 



(23) 
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The calculation is more easily performed integrating first over the variable y?, then over 9 
and finally over v. The final result is 
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The classical cold plasma limit is correctly recovered setting — )■ 0, Vp — )■ 0, 
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(27) 



yielding with u = i5 the exact expression for the growth rate 5 
1 



kVb, k ^ iOp/c, 
VbUp/ c, k:^ LUp/ c, 



(28) 



in agreement with the results from the classical cold plasma model 12 1. 

The present calculation is worth comparing to the fluid theory including a Bohm pressure 
term 0,0] and with Vp — )■ 0. In this fluid limit, the dispersion relation can be exactly solved 
and gives 



uji + cH' + ^-[[uji + c'e-^f + Aev,'uji 
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(29) 



for the unstable mode. Knowing that the fluid and the kinetic models should merge in the 
long wave length limit, we first expand the above equation for small k, 



-e + + (n^ + c^) ^ + o{k') . 

Am'' ui 
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To compare with the kinetic theory, we now expand Syy in the same limit and with 



-yy 



4m2u;6 ' (31) 
The dispersion equation (13T1) is a cubic for u"^. There are analytic formulas which can be 
applied, but in the long wavelength limit it is cheaper to solve it recursively. The result 
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FIG. 2: Typical plot of 6'^D(d), where D is the dispersion function and = + i5, for Z = 1.6 
(lower curve) and Z = 1.7 (upper curve). The intersections with the horizontal axis correspond to 
the growth rate. Parameters are /3 = 0.1 and p = 10. 

agrees with Eq. fl30|) . confirming the equivalence between quantum kinetic and quantum 
fluid models for long wavelengths. Such equivalence can be checked on Figure |3] which 
compares the kinetic with the quantum fluid with Bohm potential and the classical cold 
plasma calculations. 

In order to proceed further in the calculations, we introduce the usual dimensionless 



parameters 



20| 



c 



hwp 



(32) 

In the classical case, accounting for a thermal spread requires an additional parameter, 
namely, the temperature. Here, the thermal spread Vp is related to the density. As a result, 
the above parameters are not independent of each others and one can check that 



H 



where Oc = /he = 1/137 is the fine-structure constant. Like in the classical case 12|, the 
roots of the dispersion equation are here found with zero real part. Setting Q = + i6 and 
denoting D{6) the dispersion function yielding the dispersion equation D{6) = 0, Figure [2] 
sketches a typical plot of 6^D{6) for two values of Z. The intersection of the curves with 
the horizontal axis directly gives the growth rate. 

Figure [3] compares the kinetic growth rate obtained solving numerically the full dispersion 
equation with the classical cold and the quantum fluid with Bohm term results for the 
parameters specified in caption. The kinetic equation has been solved using Mathematica's 
"FindRoot" routine, giving the fluid growth rate ( 129|) as an initial guess. The agreement in 
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FIG. 3: (Color Online) Kinetic growth rate (blue) obtained solving the kinetic dispersion equation, 
compared to the quantum fluid result with Bohm pressure term (purple) and to the classical cold 
plasma result (yellow), in terms of the reduced wave vector Z. In (a), parameters are /3 = 0.1 and 
p = 15, and the fluid unstable range is wider than the kinetic one. In (b), parameters are f3 = 0.1 
and p = 900, and the fluid unstable range is smaller than the kinetic one. The saturation value 
for the classical cold curve is simply f3. 

the long wave length limit is clear. The classical cold plasma result saturates at the value 
P for large Z while the quantum fluid and the present kinetic results exhibit cutoffs, as the 
kinetic pressure eventually acts to prevent the pinching of small filaments. 

Note on Fig. [3b that the range of unstable modes is smaller in the fluid than in the 
kinetic case. We thus turn now to the investigation of the cutoff wave vector. 
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B. Cutoff wave vector 



The cutoff wave number Zm can be found writing = for Z = Zm and 5 = 0. 

Denoting L = lim^^o ^^D{6), one finds (only tfie numerator is shown) 

L= 4HZp[l2P\l-4p^) + Z\32 + 3H^P^p^)] (34) 



+ 3/3^ [16 - 8{H^Z^ - 8)p2 + H^Z\H^Z^ - 16)p^] In 



2-HZp 



2 + HZp 

We start searching the expression of the largest unstable wave vector Z^ in the most 
interesting limit. This is the kinetic one, with p = Vb/Vp — ?■ 0, as the opposite limit is just 
the fluid one. By developing the logarithm in Eq. fl34l) . the equation L = simplifies to 

= 384/3V^ + [3H^/3^p'^ {{H^Z'^ - 16)p'^ - 12) - 128] . (35) 

Assuming H^Z^ <^ 16 (to be checked later), we find directly 

(36) 

^/Snp + aj{3 + 4p2) 

We can now check our assumption H^Z^ ^ 16 . Replacing Z by the value of Z^, above, we 
find 

H-'Zl = < 32/3. (37) 

It thus turns out that the condition H'^Z^ <^ 16 is only weakly verified. However, Eq. f p5|) 
can be solved exactly very easily, and the solution found is numerically very close to Eq. 

m- 

Examining now the condition to expand the logarithm HZ^p <C 2, we can find the 
validity domain of Eq. 



HZmp = p\ - ^ Pfo , A 2\ '^^^ P'^ ^ ~ 38 

Equation is thus found valid in a very wide range of parameter defined by the strong 
inequality f l55]) . Note worthily, it defines various p scalings. For p <^ ac/3/87r, Z^ oc p^/^. 
Then, the denominator behaves like y^, yielding a Z^ oc p scaling until the quadratic term 
under the square root overcomes the linear one. Comparing these two terms gives a criterion 
on p almost identical to Eq. ( 138|) . 
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FIG. 4: (Color Online) Most unstable wave vector Zm in terms of p. The analytical expressions in 
the various regimes are provided on Table [H 

For p ^ 2-n /ttcP, the logarithm in Eq. f l34|) can be expanded assuming now HZ^p ^ 1- 
Performing such expansion and replacing H by its value given by Eq. (!33|) gives the following 
equation for Z^, 



8Z^ae + 9Z^a,/32 - 97r/3V(l + V) = 0. (39) 
This equation can be solved exactly. Expanding the relevant root for large p gives, 

1/4 



{Ppf/\ (40) 



Qtt 

Equations fl36ll40p eventually define three different scalings which are summarized in Table 
[H It is interesting to "unfold" the dimensionless parameters in order to explain the key 
quantity HZ p. We find, 



hUp kVb Vb hk k 

where kp is the Fermi wave number. 

TABLE I: Analytical expressions for the largest unstable wave vector in the various regimes. 



p 




Sit ^ P ^ OLcP 






/S^ 3/2 
y ac P 




{^y\ppf'' 



Figure H] displays the numerical evaluation of the cutoff wave number in terms of p and 
for 3 values of /3. The analytical expression reported in Table |T] cannot be distinguished from 
the numerical calculation within their range of validity. 
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Let us finally compare the kinetic cutoff with the fluid one. The largest unstable mode 



can be expressed exactly in the fluid model with the Bohm pressure term as 



i 



z - ^ 



^/1 + 8/hW-1 



1/2 



(42) 



Replacing H by its value in terms of p from Eq. (!33|) and expanding for small and large p's 
yield, 



3tt 
2ac 



(43) 



These extreme scalings are thus very similar to the kinetic ones explained on Table [T] in 
the limits p — )■ 0, 00. The kinetic cutoff displays three different regimes and the fluid one 
above, only two. While in the intermediate kinetic regime, the kinetic cutoff is smaller than 
the fluid one (as in Fig. [3^), Tableland Eq. (H3l) show that surprisingly, the fluid cutoff 
can be smaller than the kinetic one (as in Fig. |3b). 



C. Most unstable wave vector and largest growth rate 

The most unstable wave vector Zmax in the range [0, Zm], together with its growth rate 6m, 
are key quantities which eventually define the strength of the instability and its time scale. 
Although Fig. [1] suggests otherwise, an analytical approach based on a Taylor expansion 
near = of the dispersion equation has not been found valid here. It seems that indeed, 
the expansion up to the second order is not enough to approach the numerical calculations. 

We thus resort to a systematic exploration of the parameters phase space (/3, p) in order 
to extract scaling laws. The results are sketched on Figure [5] within the parameter range 
where numerical stability allowed to derive trustful results. 

The curves for the maximum growth rate all saturates at 6m = P for p ^ 1. In the kinetic 
regime p ^ 1, simple scalings are evidenced in terms of the parameters, and the following 
fit has been found 

6m ~ 0.72/3p2. (44) 

The largest unstable wave vector Zm, as given by Eq. ( I36l) . has been represented with 
the most unstable one Zmax on Fig. |5l We obviously find fulfilled the inequality Zm < Zmax- 
These two quantities remain locked to each other until they decouple from p* {/3) ~ 1 slightly 
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FIG. 5: (Color Online) Numerical determination of the most unstable WclV6 vector ZfYiax 

and its 

growth rate 6m, in terms of p and for various values of (3. Upper plot: Kinetic value of Zmax (plain 
curves) vs. Eq. (I36p (dashed curves). Lower plot: Kinetic value of the largest growth rate dm 
(plain curves) compared to its fluid value (dashed curves). 

varying with /3. For p <^ p*, the following equality is fulfilled with remarkable constancy, 
regardless of the value of (3 

Zmaxip < P*) ~ 0.6 Zm- (45) 

For p^ p*, Fig. |5] shows Zmax switches from a simple p to a measured p^^^ scaling. As 
expected then, the fastest growing wave number tends to infinity, but slower than the cutoff. 
A typical plot in this range is displayed in Fig. [Sti, where the growth rate quickly reaches 
its maximum ~ /3 before its progressively comes down to zero. 

Still on Fig. |5l the lower plot representing the largest growth rates shows the kinetic one 
undergoes a transition from 5m oc p to 5 m ~ /3 near p — 1. The fluid results displays a 
measured p^/^ until p ~ from where it saturates also at 5m ~ /3- If then, one wishes 
to define the limit of the fluid model through the correspondence of the maximum fluid and 
kinetic growth rates, the fluid approximation is found valid for p ^ 1. 



13 



IV. CONCLUSION AND DISCUSSION 



In conclusion, using Stratonovich's gauge invariant Wigner equation we have re-derived 
the general form of the electromagnetic dielectric tensor in a quantum plasma. The quantum 
filamentation instability was then treated using kinetic theory. The equilibrium Wigner 
function was taken as a pair of zero-temperature Fermi-Dirac distributions centered at the 
beams velocities. In this way, not only the quantum diffraction effects inherent to the kinetic 
equation ([1]) but also the fermionic character of the beams were included. The results were 
compared to the zero-temperature classical and quantum hydrodynamic equations (with 

— )■ 0), showing agreement in the long-wavelength limit. In addition, analytical expressions 
are derived for the the largest growth rate, the cutoff wave vector and the most unstable 
wave vector. 

It is worth to comment on the influence of the quantum properties against the filamen- 
tation instability. As apparent from Eq. ( 132|) . the quantum statistical effects decrease with 
the quantity p = Vt,/Vp. Also quantum diffraction effects represented by the parameter 
H = hwp/{mV^) decrease with p, as follows from Eq. f l5^ for a fixed beam velocity. On 
the other hand, from Figs. dH E]) we see that the cuttof wavenumber, the largest unstable 
wave vector and its growth rate increases with p. Hence kinetic theory shows that quantum 
mechanics has a stabilizing role to the filamentation instability. The present theory can be 
applied to very dense counter-streaming plasmas obeying Fermi-Dirac statistics and where 
quantum diffraction can be significant. 
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